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I. INTRODUCTION 


The contract was initially oriented toward simulating the radar 
parameters of the ASAR (Advanced Synthetic Aperture Radar) system under 
development jointly by NASA/JPL and NASA/JSC. The former had been 
tasked to develop a system using a technique known as "burst mode" 
with wavelength diversity in pulse bursts of transmissions. NASA/ JSC 
was responsible for developing the antenna system required for handling 
the wide band of frequencies Involved. The ARL:UT supporting effort 
was to simulate the system design and particularly the "burst mode", 
to provide insight into the operation of the algorithms and anticipated 
results. The simulation developed by ARL:UT for an orbiting SAR 
platform had to be modified to support the aircraft based ASAR. 

After three months of effort toward these objectives, a decision 
was made to stop further work on the ASAR system and Instead concentrate 
on certain key SAR technology areas. Among these are multi,.olarlzation 
antennas, multifrequency radar, multibeam squint mode SAR (wide swath), 
calibration techniques, and burst mode implementation. Of these, 

ARL:UT addressed the first three in varying degree. 
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II. ASAR SIMULATION EFFORTS 


A. Discussion 


The ARL:UT orbital SAR simulation (OSS) had to be modified to 
simulate an aircraft flying in earth's atmosphere at a relatively low 
velocity and altitude compared to the orbiting platform. In addition, the 
radar parameters had to be obtained from the systr i designers at NASA/JPL 
and the antenna parameters from the project personnel at NASA/JSC. The 
system was to have exhibited two unique design features: extremely 

broad wavelength coverage, and data processing of synthetic data arrays 
formed by bursts of pulses at several different radar frequencies. 

The purposes of the simulation were to investigate and demonstrate 
the validity of the system design, in particular the so-called "burst 
mode” operation. 

Two sets of problems were encountered in this effort. The first 
had to do with adapting the existing simulation, of a space-based 
orbiting radar, to the aircraft-based design. This problem set was 
rather easily solved, but attempts to define the design parameters of 
the radar system were relatively unsuccessful initially, and two or 
three months elapsed before a preliminary set of parameters for use in 
the simulation was agreed upon and obtained. 

B. Dynamics Simulation 

In order to utilize the existing routines in the simulation, 
the proposed aircraft platform (a CV 990) was considered like an 
orbiting platform in a low earth gravity environment, so that 
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it could maintain orbit at Its usual altitude and velocity. The 
magnitude of the velocity vector Is 

V - a(l-e^cos^E)^^^ * E , (D 

8 

with the altitude specified by the orbit radius from earth's center 


r^ “ a(l - e cosE) 


( 2 ) 


The gravitational constant of earth u has been measured at 
398,601 km /sec . In Eq. (2) e Is the eccentricity of the orbit, and E 
Is the eccentric anomaly. If e Is zero, the orbit Is circular. The 
parameter a Is the semlmajor axis of the orbit ellipse, or the radius 
of the orbit with e«0. 


By selecting the proper values for the six Kepler orbital 
elements, a, e, 1, 0, u, and T^, the orbit Is specified. In Eq. (1), 

E Is the time derivative of the eccentric anomaly: 

E - e sinE - v(t-T ) (3) 

P 

•j— (1 - e cosE) = V . (4) 

ut 

Here v Is the mean angular velocity of the platform 

V « J \i . (5) 

P 

In order to derive a usable value of y we can simplify the 

P 

operation by assuming a spherical earth and a circular orbit for the 
aircraft. We find that 

r*a , V«a»<E , E*v ; (6) 

s s 
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therefore 


v= - . 

Assuming that the CV 990 will fly at 400 kt TAS (205 m/sec) at an 
altitude of 12,000 m (39,372 ft), we can solve for an Imaginary value 
of which will provide the proper angular rate. 

y * avf (8) 

p s 

- (0.205)^ (1246378.167) 

“ 268.2 km^/sec^ , 


compared with the measured value for earth cf 

y - 398,601 km^/sec^ 

P 

All other planet earth parameters remain the same. 

For the orbit specification, all parameters used for SEASAT, for 
example, can remain the same except the semimajor axis a, which becomes 
6390.167 km (1246378.167) and the orbit rotational rate which becomes 


. 1 - ( 0»205 \ 

2ir r 2 it V6390.167/ 
s 


18.4 >« lO”^ deg/sec 


C. Ra<iar Simulation 


The proposed ASAR design was to have used linear FM (chirp) pulse 
compression, whereas the OSS Is designed for binary phase coded pulse 
compression, which Is particularly adapted to digital processing. 
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Some effort was spent in analyzing the pulse compression routine 
requirements to provide an analogous digital signal sampled at the 
Nyqulst rate. 

Following a specification review at JSC, a first attempt was 
made to assign radar specification values for the simulation. Table I 
contains most of these values. 


The linear FM chirp modulation proposed for the ASAR has the 
following parameters. 

f - 10® HZ ; f - 1.02 X 10® Hz . 4f " 10^^ Hz/sec 
o 1 at 

T - 20 X 10“^ sec 


The equation for the phase as a function of time is 
att “0 , 

o 

at t » n , 4 » » <► + f (n) + x ^ (n)^ 

n o o z at 

at t » n -1 , “ '•'o "2 ^ 

If At - n - (n-1) - 1 . A* - ^ ^ ^ (2n-l) . 

From the sampling theorem, sampling period " 1/2 period of 
highest frequency, 

let At “ - 57 “ (n-(n-l)] - - 5 ^ 5 0.49 nsec , (10) 

^'1 ^'1 2.04 X 10 ^ Hz 


then 

♦n - +0 + 2 ^ 


( 12 ) 



TABLE I 


Radar Specifications Radar ID; ASAR 1 


L S C X 


X 

WL 

Radar wavelength, m 0.23513 0.093685 0.062457 0.0*^ 

TB 

TBW 

Range time-bandwidth product 

400 



S/N 

Received signal-to-noise ratio, dB 

7.0 


N 

P 

LRC 

Sample length of range correlation 

40.800 


% 

LPI 

Sample length across phase Interval 

1 


Pr 

RESR 

Ground range resolution, m 

25 y/l 


Pa 

RESA 

Azimuth resolution, m 

2 


s 

r 

SRR 

Range sampling ratio 

1 



KPC 

Phase code, 1> -10“\[o.98 + (2n-l) 
n 

0.24 X lO' 

•3] r.d 


KPCS 

Phase code sequence, BCD 

1 
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a 

SRA 

Azimuth sampling ratio 

1 



KRFTNR 

Name of range Impulse response function, BCD: 

cosine^ 


KWTFTNA 

Name of aperture weighting function 

f BCD 

Taylor 


SHADFAC 

Aperture weight factor 


-30 dB 

POR 

PCOFFR 

Patch-to-patch offset, range, m 


0 

POA 

PCOFFA 

Patch-to-patch offset, azimuth, m 


0 

SW 

SW 

Swath width, km 


20 


NA 

Number of patches 


1 

LAT 

P 

STLAT 

Map start latitude, rad 


2n/10 


TCR 

Planetocentrlc distance to terrain ( 

center, km 


t 



0.63675192370 x 10’ 

LT^ 

TCLAT 

Terrain center latitude, rad 


2n/10 

LNt 

TCLONG 

Terrain center longitude, rad 


ir 
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1 df 29 9 

■ Vl * '^*=^0 2 dt ^ + (2n-l)Atn 


- Atf^ + 2 dt <2n-l)At^ cycles 
■ 2irf^At + IT ~ (2n-l)At^ rad 


and 


- 2ir(10^)(o.A9xlo"®) + tt . (2n-l)(o.49xlo"® 

' 20 X 10“° ' 

= 10"^n 1^0.98 + (2u-l)(o.24 10”^)jrad 

- 0.18 [0.98 + (2n-l) (0.00024)] deg 

N - TT “ (2.04x10®) (20xi0“^) - 40,800 samples, 
p At 


0. Antenna Simulation 


(13) 


The antenna design was the iresponslblllty of JSC. It posed a 
special problem because a single antenna for the 10 GHz bandwidth 
(2-12 GHz) would be very difficult to design and build. Secondly, Its 
design largely depended on the JPL radar design, which was not finalized. 
Table II gives the parameters that were selected for use In the absence 
of any real antenna characteristics. 

E. Terrain Specifications 

The terrain to be used Initially was to be a model terrain used 
In previous work, called SINGLESCA'.. It contains two discrete scatterers 
and a single homogeneous field. Table III describes the parameters for 
this model. 


8 



TABLE II 


* 


KANT 

ANTENNA SPECIFICATIONS 
Antenna identification, BCD 

UNITPAT 

* 

"b 

SNAC 

Boresight nadir at oap start, rad 

20*,ir/9 

* 


SQUINT 

Boresight squint at map start, rad 

90*,x/2 

* 

az 

ABW 

Anteana azimuth angular coverage, rad 

w 

* 

el 

EBW 

Antenna elevation angular coverage, rad 

IT 

* 

PHCPB(3) 

Phase center position, body axis xyz, m 

0,0,0 

s 


ABRPYO) 

Antenna attitude, body axis rpy, rad 

18*4 

* 


PDRPY(3) 

Platform attitude rates, local orbital 
rpy, rad/sec 

0,0,0 
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TABLE III 


TERRAIN SPECIFICATIONS 


KSNO 

Terrain identification, BCD 

SINGLESCAT 

NDISC 

Nuaber of discretes in terrain aodel 

2 

NFLDS 

Ninber of fields in terrain wodel 

] 

NSCATS 

Total nuaber of scatterers in terrain nodel 

102 

DffiAN 

Mean echo strength 

0.1 

DISTAX 

x-axis coverage, ka 

10 

DISTY 

y-axis coverage, ka 

10 
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A seetlng was held at NASA/JPL in the third nonth of the contract, 
attended by representatives fron ARL:UT and Environnental Research 
Institute of Michigan (ERIM), to evaluate the design objiectives and 
status of the ASAR. The project had schedule and funding difficulties 
and, subsequently, NASA decided to cancel the entire project. Unfortunately, 
a substantial portion of the contract resources at ARL:UT had already been 
expended on the effort, and a sioulation had not been run. 

The project was reoriented with specific eaphasis on three areas: 
antenna polarization effects, wavelength diversity effects, and iapleaen- 
tatlon of a so-called "burst aode" wavelength diversity concept on a 
patch-by-patch basis. 
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III. OTHER SIMULATION INVESTIGATIONS 

Figure 1 illustrates the architecture of the OSS siaulatlon. Each 
interaediate step results in data on tape or disk bulk storage, which 
peralts increaental coapletion of the entire simulated radar laage. For 
the polarization and wavelength diversity ezperiaents an existing terrain 
aodel was utilized, illustrated in Fig. 2. This figure Illustrates the 
scene content of the test terrain as well as laages derived froa it. 

A. Antenna Polarization Effect 


A coaplete siaulatlon would Include the polarization sensitivity of 
individual scatterers and would incorporate the randoaly polarized 
backscattering froa the individual scatterers in fields of scatterers, 
particularly hoaogeneous fields. The general situation is that the 
backscattered wave is partially polarized, part of it coapletely polar- 
ized, and part of it unpolarlzed. Suppose we designate the degree of 
polarization as the ratio of the polarized power to the total power 
(following Kraus), or: 


j „ polarized power ^ 
total power 


(15) 


and denoting the Stokes paraneters as s^, s^, 82, s^, 




2 ^ 2 ^ 2 

®1 ®2 ®3 


0 < d < 1 


(16) 
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These are normalized Stokes parameters > defined as follows , with 
reference to Fig* 3: 

Ex » £ cosc sinu)t cost - E sine cosa>t sint (17 

o o 

Ey » E cose sinu)t sint ^ E sine cosoit cost • (Ifi 

o o 

If we now substitute the time dependent values of Ex and Ey we 


E« sin((i)t-6«) * E (cose sina>t cost - sine cos(ot sinT) (19) 

1 1 o 

E^ 8in(u)t-5^) « E (cose sin:^t sinr + sine cosojt cost) • (20) 

2 2 o 

If these are expanded, and sincot terms and cosmt terms set equal (which 
eliminates the tic;e dependence): 


■ "o(‘ 


2 2.^2 ^21 
cos e cos T + sin e sin t. 


( 2 2 2 2 
cos e sin t + sin e cos^rj 

2 — 
The flux density in W/m is the Poyntlng vector S: 


s « s + s - 

' ' X y 


2 2 2 
^1 + ^2 


where is the free space impedance, 377 Q per square measure, is the 

Poyntlng vector for the E wave component polarized to x, and S , for the 

o y 

one polarized to y. 


*^1 ”” 2 2 2 2 
S « • |s|(co8 e cos t + sin e sin t) 

o 


o 2 2 , 2.^2 2 . 

S ■ » |S|(cos e sin t + sin € cos t) 

The Stokes parameters are now defined as 


1 

^ 

o S 


cos2c cos 2t 
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FIGURE 3 

POLARIZATION ELLIPSE 
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where <> denotes time averaged value. Then 


^2 * 'sz ^2 <cos2e sin2x> 


S3 » — <E^ cos5> = <sln2e> » 

2 2 ^ 

1 > + S3 

These can be written in matrix form 


■ = s. 


i - 0,1, 2 , 3 


0 d cos2e cos2t 

= S + S ( 27 ) 

0 d cos2e sln2t 

0 _d sln2e 

for a partially polarized wave, where the first term of Eq. ( 27 ) is the 

polarized power density and the second Is the unpolarized power density. 

Tlie antenna polarization response may also be expressed as a matrix, 
with the effective aperture: 


“1 2 

A [a . ] “ A m 

el e 


Thus the power out of the antenna Is given by 


W - 2 S (sj 


2 ^ ®i 
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The simulation for Che APQ-102 uses a linear polarized antenna, 


with a gain given by 


4 IT A 



( 30 ) 


and the magnitude of the gain function can be entered into Che simulation 
via program ANTENA. Substituting the matrix representation for the A^ 
we have 


0[.jl . 


X e 


^2 


(31) 


Combining the polarization matrix for a scatterer, we have 




4TL ioi r 

x2 2 ^ 


Vl®2®3^ 




(32) 


where the substitution of radar cross-section a is made for the power 
density S, assuming that the backscattered power density la proportional 
to the cross section. 


To implement this function in program ANTENA requires that the a^ 
be specified, as well as the s^. The former may reduce to 
82 * 0 , 8^*0, or it may depend on platform motion. 

The specification of [s] on the other hand would be arbitrary or 
would depend on measured data for particular types of radar targets or 
clutter. 
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As a first step, it /as decided to investigate the effects of 
crosspolarized antennas for transmit and receive, and assume no polariza- 
tion sensitivity for a. Although the crosspolarized response of the 
receive antenna has been measured, it was not available in the azimuth 
plane, only in elevation. Based on the lobe structure and boresight 
gain of the elevation pattern, a crosspolarized azimuth response pattern 
was constructed and used to specify ANTENA. Figure 4 illustrates the 
pattern used. 

The values for the patterns, digitized every 0.33® to a coverage of 

124°, were entered using an HP9810A with digitizer. The SEASAT radar 

parameters were used and six scatterers were entered into the terrain 

model for test purposes. Five of these were given a cross section of 
2 

1 m at the center and corners of the 1.4 x 1.4 km mapped area; the other 

2 

was 26 dB greater at 20 m . The imulation was exercised and the SARCON 
data were printed out — see the Appendix. Figure 5 shows the displayed 
results for the co-polarized antenna response. Post-processing was 
effected at 2 dB/gray shade, with the peak filter magnitude falling in 
gray shade 16. Figure 5(a) has gray shade 16 turned on and Fig. 5(b) 
has it off. To expedite processing no range pulse compression was 
used; hence range sidelobes are not present. Figure 6 is for the same 
scatterer array but uses a unit antenna pattern, i.e., it is isotropic. 
There is no appreciable difference since the entire imaged area falls 
well within the beam of the APQ-102 1.3° 3 dB beamwidth for the orbiter 
altitude of 870 km. 

To demonstrate antenna response, a realistic terrain scene without 
a prohibitively large scattering area was required and it was necessary 
to place scatterers at ambiguous (with respect to pulse repetition 
frequency (PRF)) response points. The total azimuth coverage imaged 
was only 1.4 km; if the terrain model covered only this area, then from 
an altitude of 870 km along 4.6 km of the orbit (the synthetic array 
length), only ±0.25° of the antenna patterns would be used. To create 
an effective but simple model of terrain, single scatterers were placed 
at ambiguous points as Illustrated in Fig. 7. Only the single scatterer 
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FIGURE 4 

APQ-102 ANTENNA SIMULATION PATTERNS 
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(a) Six-scatter OSS simulation* 
a 1 1 gray shades , br ifhtoess on 

display set for mlniiwon Doppler 
si delobes on 100 x scatterer 


(b) Gray shade 16 turned off, 

brightness increased - to incl ude 
all Doppler filter sidelobes 


FIGURE 5 

OSS SIMULATION OF SIX-SCATTERER ARRAY USING SEASAT SYSTEM 
WITH APQ-102 ANTENNA PATTERN IN AZIMUTH, 

NO PULSE COMPRESSION 


AE-81-133 


ORIGINAL PAGE 

BLACK AND WHfTE PHOTOGRAPH 



(a) Six-scatter OSS simulation, 
all gray shades, brightness on 
display set for mini mum Doppler 
si delobes on 100 x scatterer 



(b) Gray shade 16 tur.ied off, 

brightness Increased to include 
all Ooppler filter si delobes 


FIGURE 6 

OSS SIMULATION OF SO-SCATTERER ARRAY USING SEASAT SYSTEM 
WITH UNIT ANTENNA PA'TTERN IN AZIMUTH, 

NO PULSE CWPRESSIW 


AE-81-134 
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PATCH CENTER 

27 Discrete Scatte'-ers in Model 


FIGURE 7 

AMBIGUITY TERRAIN 
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at the patch center was imaged; however, there were Doppler contributions 
from all the sidelobe scat ter ers. Two configurations of scatterers 
were used; the second one offset the ambiguous point scatterers in 
graduated steps 10-30 m from the ambiguity locations. 

Figure 8 illustrates the results* Each photograph contains six 
images of the patch center; each of these is 1*4 x 0*4 km, with 25 m 
resolution* From «:he top, the first three strips are of the unoffset 
terrain model imaged with the isotropic, co-polarlaed, and cro6spo?ariaed 
patterns. All three patterns were normalized to 1*0 at the peak response* 
The 1 . it three from the top on the left are images of the offset scatter- 
ers for co-polarized, crosspolarized, and isotropic patterns* The top 
gray shade of the display is set down 20 dB from the peak response of 
the center pixel to show the weak sidelobe structure* 

The right photo has the same order of images, but the top gray 
shade Is 35 dB below the peak response, thus defining more clearly 
sidelobe response of the system configurations* 

Analysis of these results shows a significant reduction in 
ambiguities between images with the isotropic antenna pattern and the 
APQ-102 patterns. The co-polarized pattern image shows no ambiguities 
and has lower sidelobes than the crosspolarized image* The cross- 
polarized pattern generates an ambiguity due to its high sidelobes « One 
must Iz careful in trying to interpret the power level of the ambigui- 
ties, because the energy returned outside the main beam is from only a 
few scatterers and may be lover than the energy returned from an actual 
continuous terrain area. 

B. Frequency Diversity Effects 

It was decided to investigate the image variations with wavelength 
on the scene of Fig. 2, which h.is been imaged at the SEASAT wavelength 
of 23.5 cm. Three other waveit ^ths were selected: 1.8, 3.125, and 

10.34 cm. 
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To properly implement the simulation the variation of the cross 
section with wavelength should be included. Other wavelength dependent 
factors are the antenna response patterns, transmitted power, receiving 
system noise and. Indeed, every factor in the radar range equation either 
implicitly or directly. In the interest of economy and reduced com- 
plexity, however, only the effects on the synthetic array processing 
were taken into account. The results are not very dramatic, as would 
be the case if, for example, the cross-section dependency were Included, 
This could involve a simple adjustment since, generally speaking, radar 
cross section varies proportionally to the square of the wavelength, A 
useful set of relations might be: 

2 2 
* 1 m =0 dB m 

0j^(lO.34)^ 

^ = 1. 936 X lo"^ o, = -7 dB 
^ (23.5)^ ^ 

Ov = ^ = 1.768 X 10 ^ a, = -17,5 dB m 

(23.5) ^ ^ 

" T “ 5-867 X 10"^ o, = -24 dB 

(23.5) ^ L 

This is in general true only for perfectly conducting bodies, and 
is not the general case for clutter targets. For the first simulation 
effort, only the wavelength parameter was varied, and the resolution 
cell dimensions were forced to the same values by varying the number of 
pulses processed or, equivalently, the array length. 

Figure 9 presents the four images of the test pattern, photographed 
with two different exposures. The most obvious effect of the wavelength 
changes is in the coherent speckle pattern " >r the hotaogeneous fields. 
The data were analyzed statistically, and only very slight changes In 
the average values were observed. Table IV provides these data for 
comparison. The simulations were run for the SEASAT radar and the 
isotropic antenna pattern. The ARLiUT high resolution display system 
was used to display and photograph the results (Fig, 10), 
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TABLE IV 

PARTIAL RESULTS OF FREQUENCY DIVERSITY STUDY 
STATISTICAL DATA 



Wavelength 

Filter Magnitude 

Humi 



liHi 

: Mexlnum 

336.631A7 

337.97644 

334. <7785 

342.83377 

1 

Minimum 

-12 

6.86>‘10 

-14 

3.69x10 

-13 

5.13x10 

8.95x10"^^ 

Mean 

0.18181 

0.17562 

0.17848 

0.180995 

Sigma 

4.857 

4.859 

4.849 

4.959 

i 


SYNTHETIC ARRAY PARAMETERS 


Array No . 1 

43,262 pts 

43,039 pts 

42,816 pts 

42,816 pts 

Array Length, m 

4602 

2026 

612 

352.5 

Formation Time 
(sec) 

0.617 

0.271 

0.082 

0.047 

No. of Pulses 

926 

916 

900 

896 

PRF. Hz 

1,500 

3,370 

10,970 

18,948 

Platform Velocity 
(km/sec) 

7.4577 

7.4577 

7.4577 

7.4577 


The PRF was increased to obtain approximately the same number of 
samples (pulses) in array lengths progressively shorter so that resolu- 
tion would remain the same for the aximuth dimension as wavelength 
decreased. Had this not been done, the azimuth resolution would have 
drastically increased as the ratio of array length to wavelength Increased. 
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IMAGE CONTOURS FROM 
SELECTED GREY SHADES 


FIGURE 10 

ARL:UT HIGH RESOLUTION DISPLAY SYSTEM 


2838-14A 


W: 






IV. TECHNICAL EFFORT SUMMARY 


Unfortunately, the dilution of the resources allocated to the 
simulation efforts by the ASAR project prevented completion of in-depth 
studies on any of the subjects undertaken. Initial analysis was com- 
menced on a wide swath technique using squinted multiple beams. A major 
challenge will be the formation of the image from the output of multiple 
beams tracking the clutter along the velocity vector. Figure 11 illus- 
trates the concept, which will be undertaken in a follow-on effort. 

In summary, the ASAR work provided a method for adapting the OSS 
simulation to an aircraft platform, and to a linear FM (chirped) pulse 
radar. The analysis indicates the approach to be taken to fully char- 
acterize the effects of polarization of the backscattered energy and 
crosspolarization of the antenna system versus co->polarization of trans- 
mit and receive antennas. 

Finally, some effects of frequency diversity in the generation of 
coherent speckle from homogeneous fields were synthesized with the simu- 
lation* These efforts point the way to more detailed and fruitful 
investigations in the future, using the OSS simulation programs. 
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RECEIVE BEAM 
COVERAGE 


FIGURE 11 

THE GEOMETRY OF THE MULTIPLE BEAM SAR 


ARLiUT 

AE-81-137 

CRG-GA 

9-16-81 




L-BAND WAVELENGTH DATA 
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